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Abstrac t 

The results of an analytical and experimental study 
of the nonlinear response of thin, unstiffened, aluminum 
cylindrical shells with a long longitudinal crack are pre- 
sented. The shells are analyzed with a nonlinear shell 
analysis code that accurately accounts for global and lo- 
cal structural response phenomena. Results are presented 
for internal pressure and for axial compression loads. 
The effect of initial crack length on the initiation of stable 
tearing and unstable crack growth in typical shells sub- 
jected to internal pressure loads is predicted using geo- 
metrically nonlinear elastic-plastic finite element 
analyses. The results of these analyses and of the exper- 
iments indicate that the pressure required to initiate stable 
tearing and unstable tearing in a shell subjected to inter- 
nal pressure loads decreases as the crack length increases. 
The effects of crack length on the prebuckling, buckling 
and postbuckling responses of typical shells subjected to 
axial compression are also described. For this loading 
condition, the crack length is held constant. The results 
of the analyses illustrate the influence of crack length on 
shell buckling instabilities. The experimental and analyt- 
ical results correlate well. 

fott odustion 

Transport fuselage shell structures are designed to 
support combinations of internal pressure and mechani- 
cal loads which can cause the structure to have a geomet- 
rically nonlinear structural response. These shell 
structures are required to have adequate structural integ- 
rity so that they do not fail if cracks develop during the 
service life of the aircraft. The structural response of a 
transport fuselage shell structure with a crack is influ- 
enced by the local stress and displacement gradients near 
the crack and by the internal load distribution in the shell. 
Local fuselage out-of-plane skin displacements near a 
crack can be large compared to the fuselage skin thick- 
ness, and these displacements can couple with the inter- 
nal stress resultants in the shell to amplify the magnitudes 
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of the local stresses and displacements near the crack. In 
addition, the stiffness and internal load distributions in a 
shell with a crack will change as the crack grows and 
when the skin buckles. These changes will affect the lo- 
cal stress and displacement gradients near the crack. 
This nonlinear response must be understood and accu- 
rately predicted in order to determine the structural integ- 
rity and residual strength of a fuselage structure with 
damage. Recent studies (e.g., Refs. 1-5) have shown that 
the structural response and structural integrity of a shell 
with a crack can be studied analytically with a nonlinear 
structural analysis procedure that can model crack 
growth in the shell. 

Typical nonlinear analysis results presented in 
Refs. 3 and 4 indicate that different combinations of ap- 
plied loads can cause different responses of a stiffened 
shell with a long crack. The magnitudes of the mechan- 
ical loads used in the studies described in Refs. 3 and 4 
are representative of loads that do not buckle the skin of 
the fuselage. Fuselage shells are usually designed to al- 
low the fuselage skin to buckle above a specified design 
load that is less than the design limit load for the shell. 
During the design of the fuselage, it is assumed that the 
design limit load can occur anytime during the service 
life of the aircraft. As a result, a long crack could exist 
in the fuselage shell after a considerable amount of flight 
service, and loading conditions could occur that cause 
the shell with the long crack to buckle. The results of a 
preliminary analytical study 5 of the effects of long 
cracks on the nonlinear response of unstiffened alumi- 
num shells indicates that the behavior of a shell can be 
influenced significantly by the initial length of the 
crack . 

The present paper describes the results of an analyt- 
ical and experimental study of the effects of internal 
pressure and axial compression loads on the responses of 
thin, unstiffened, laboratory scale, aluminum cylindrical 
shells with a long longitudinal crack. Two wall thick- 
nesses are considered in the analytical study to determine 
the effect of wall thickness on the shell response. The 
predicted effect of the initial crack length on the initia- 
tion of stable tearing and on unstable crack growth is dis- 
cussed for shells subjected to internal pressure loads. In 
addition, predicted prebuckling, buckling, and initial 
postbuckling results are presented and compared for cy- 
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lindrical shells subjected to axial compression loads. 
The results presented illustrate the influence of the load- 
ing condition and crack length on shell crack-growth in- 
stabilities, and on shell buckling instabilities. 

Shell Model and Analysis Procedure 

Shell 

The geometry of the shells analyzed in this study is 
defined in figure la. The shells have a 9.0-inch radius, 
R , a 0.040- or 0.020-inch-thick wall, t , and a 36.0-inch 
length, L . A longitudinal crack is located at 9 = 0° 
and at shell midlength. The initial crack length, a , rang- 
es from 1.0 to 4.0 inches. The shells are typical labora- 
tory-scale cylindrical shells and are made of 2024-T3 
aluminum alloy. The Young’s modulus, E, for the alu- 
minum alloy is equal to 10 msi and Poisson’s ratio is 
equal to 0.3. The yield stress for the material is 50 ksi 
and the ultimate stress is 72 ksi. 

Non lin ear An alysi s Pr ogsdure 

The shell responses were studied numerically using 
the STAGS (STructural Analysis of General Shells) non- 
linear shell analysis code. 6 STAGS is a finite element 
code designed for the static and dynamic analysis of gen- 
eral shells and includes the effects of geometric and ma- 
terial nonlinearities in the analysis. The code uses both 
the modified and full Newton methods for its nonlinear 
solution algorithms, and accounts for large rotations in a 
shell by using a co-rotational algorithm at the element 
level. STAGS has static and transient analysis capabili- 
ties that can be used to predict local instabilities and 
modal interactions that occur due to destabilizing me- 
chanical loads, such as an applied compression or shear 
load. The Riks pseudo arc-length path following 
method 7,8 is used to continue a solution past the limit 
points of a nonlinear response. With this strategy, the in- 
crementally applied loading parameter is replaced by an 
arc-length along the solution path which is then used as 
an independent loading parameter. The arc-length incre- 
ments are automatically adjusted by the program as a 
function of the solution behavior. 

STAGS can also perform crack-propagation analy- 
ses, and can represent the effects of crack growth on non- 
linear shell response. A nodal release method and a load- 
relaxation technique are used to extend a crack while the 
shell is in a nonlinear equilibrium state. 2 The condition 
for crack extension is based upon a fracture criterion. 
When a crack is to be extended, the forces necessary to 
hold the current crack tip nodes together are calculated. 
The crack is extended by releasing the nodal compatibil- 
ity condition at the crack tip, applying the equivalent 
crack-tip forces, and then releasing these forces to estab- 
lish a new equilibrium state, which corresponds to the 
longer crack. The changes in the stiffness matrix and the 
internal load distribution that occur during crack growth 
are accounted for in the analysis, and the nonlinear cou- 


pling between internal forces and in-plane and out-of- 
plane displacement gradients that occurs in a shell are 
properly represented. Output from STAGS, associated 
with a crack, includes the strain-energy-release rate in an 
elastic analysis, and the crack-tip opening angle (CTOA) 
in an elastic-plastic analysis. These quantities can 
then be used as part of a fracture criterion in an elastic 
analysis or an elastic-plastic analysis to predict the resid- 
ual strength and stable crack growth behavior of a dam- 
aged shell. 

The axial compression loads and internal pressure 
loads were modeled in STAGS using two separate load 
sets. For the case of axial compression only, the com- 
pression load was applied to the ends of the shell by 
specifying a uniform end displacement. For the internal 
pressure load cases, internal pressure was simulated by 
applying a uniform lateral pressure to the shell wall and 
applying an axial tensile force to the ends of the shell 
with multi-point constraints to enforce a uniform end dis- 
placement. The loading conditions were applied to the 
finite element models shown in figure 1 . Both the stan- 
dard 410 element and the 5 10 transition element from the 
STAGS element library were used in the models. Radial 
displacement constraints were imposed on the ends of 
the shells to insure that the ends of the shells remain cir- 
cular during the load application. 

For the internal pressure load case, advantage was 
taken of the symmetry of the problem, and only a quarter 
of the shell was modeled, as shown in figure lb. Sym- 
metry conditions were applied along the edges 0 = 0° 
and 0 = -180° and along the edge x = 18 inches. 
Cracks with initial half-lengths of 0.5 to 2.0 inches were 
defined in the model along 0 = 0° and extended from 
the symmetry boundary at the edge, x = 18 inches. 
Growth of the cracks was studied for a quasi-static in- 
crease in the magnitude of the internal pressure loading. 
An elastic-plastic analysis, with a critical crack-tip open- 
ing angle (CTOA) fracture criterion, 9,10 was used to sim- 
ulate the stable crack growth in the shell. In addition, 
straight cracks were assumed, so that crack extension 
was self-similar. The White-Besseling theory, with a 
piece- wise linear stress-strain curve, 11 was used to rep- 
resent the material nonlinearity. A high level of mesh re- 
finement was required around the crack tip to predict 
accurately the yielding at the crack tip, and crack exten- 
sion. Crack extension was assumed to occur when the 
crack-tip opening angle was equal to 5.36°. 

The entire cylindrical shell was modeled for the ax- 
ial compression load case. The finite element model for 
this load cases is shown in figure lc. Straight cracks with 
lengths of 1.0 to 4.0 inches were defined in the model 
along 0 = 0° and at the shell midlength. Plasticity ef- 
fects were included for some of the analyses for this load 
case. The primary objective of the study for the axial 
compression case was to continue the analysis beyond 
the critical buckling state, and to focus on the postbuck- 
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ling response of the shell with a crack. Thus, the interac- 
tion between local buckling and crack extension was not 
specifically addressed; that is, the crack lengths were 
held constant throughout the analyses. 

The prebuckling, buckling and postbuckling re- 
sponses of the shells for the axial compression loading 
condition were determined using the following analysis 
procedure. The prebuckling responses were determined 
using the geometrically nonlinear quasi-static analysis 
capability in STAGS. The Riks pseudo arc-length path 
following method, described above, was used to com- 
pute the initial shell response just after buckling. The 
initial, unstable, postbuckling response of the shell was 
predicted using the nonlinear transient analysis option of 
the code. The transient analysis was initiated from an 
unstable equilibrium state just beyond the buckling point 
by incrementing the end displacement. The transient 
analysis was continued until the transient response 
damped out or decayed. Once the transient analysis con- 
verged to a steady-state solution, the load relaxation op- 
tion of the code was used to establish a stable equilibrium 
state. The subsequent stable postbuckling response of 
the shell was computed using the standard nonlinear, 
static analysis option. 

Fracture Criterion 

As mentioned in the previous section, crack exten- 
sion for the pressure loading case was simulated using 
the CTOA fracture criterion. The CTOA fracture criteri- 
on uses the crack opening angle, measured at a fixed dis- 
tance from the crack tip, as the fracture parameter. The 
criterion assumes that stable crack growth will occur 
when the angle reaches a critical value. The critical an- 
gle for a particular material and thickness can be ob- 
tained by numerically simulating the fracture behavior of 
a laboratory specimen, using an elastic-plastic finite ele- 
ment analysis, and determining the angle that best de- 
scribes the experimentally observed fracture behavior. 
The determined angle can then be used to predict the 
fracture behavior of different structural configurations. 
In the numerical simulation, a characteristic distance be- 
hind the crack tip is defined for evaluating the angle. The 
crack opening displacements behind the crack are moni- 
tored and the crack is advanced by one element length 
when the opening angle evaluated at the characteristic 
distance reaches the critical value. In the present study, 
the characteristic distance was equal to 0.04 inches (min- 
imum element size along the line of crack extension). 
Newman 9 has shown this value to be adequate for ana- 
lyzing stable tearing in a wide variety of materials. The 
critical angle of 5.36° used in this study, is a preliminary 
value that was determined by matching STAGS predic- 
tions for the fracture behavior of a shell with the geome- 
try described above and with a 4.0-inch-long 
longitudinal crack, with the experimentally observed be- 
havior of an aluminum shell of the same geometry and 


subjected to internal pressure. 

Experiments 

Test S pgg im sns 

The specimens tested in this investigation were fab- 
ricated from 0.040-inch-thick 2024-T3 bare aluminum 
alloy sheet. All specimens were 39 inches long and 18 
inches in diameter. Each specimen was fabricated with 
a 1.5-inch-wide double lap splice with 0.040-inch-thick 
splice plates and a single row of 0.1875-inch-diameter 
aluminum rivets on each side of the splice. The rivets 
were spaced 1.5 inches apart along the length of the spec- 
imens and were located five rivet diameters from the 
edges of the splice plates. Longitudinal cracks were ma- 
chined in the specimens at the mid length of the speci- 
mens, and the cracks were diametrically opposed to the 
splice in the specimen. The cracks were 0.010 inches 
wide and either 2.0, 3.0 or 4.0 inches long. Specimens 
subjected to internal pressure loads had the crack tips 
sharpened with a sharp knife. 

Apparatus and_ Test s 

Both ends of each specimen were potted in an alu- 
minum-filled epoxy resin to assure proper load introduc- 
tion during the tests. The potting material extended 1.5 
inches along the length of the specimens at each end, and 
the test section of the specimens was 36 inches long. 
Specimens subjected to compression loads were ma- 
chined flat and parallel, and tested in a 300-kip hydraulic 
testing machine. Specimens subjected to internal pres- 
sure loads were attached to end fixtures, and slowly pres- 
surized to failure with nitrogen gas. One end of the 
pressurized specimens was unconstrained in the axial di- 
rection to allow free expansion of the specimen. All 
specimens were instrumented with electrical resistance 
strain gages and direct current differential transducers 
were used to measure displacements. Crack wire gages 
were applied at a crack tip of the specimens subjected to 
internal pressure loads to record the stable tearing of the 
crack. A moire interferometry procedure was used with 
the specimens subjected to axial compression loads to 
observe the deformation patterns before and after buck- 
ling had occurred. Three non-collinear direct current dif- 
ferential transducers were used to measure the 
displacements of the loading platen for the compression 
loaded specimens. All data were recorded with a data ac- 
quisition system, and the response of all specimens was 
recorded on video tape. 

Results and Discussion 

The nonlinear analysis and test results for thin un- 
stiffened aluminum cylindrical shells with a longitudinal 
crack are presented in this section. Results have been 
generated for two loading conditions: internal pressure 
only, and axial compression only. Results for these load- 
ing conditions are presented for shells with a longitudinal 
crack at shell midlength and with initial crack lengths of 
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hO, 2.0, 3.0 and 4.0 inches. The maximum value of the 
applied internal pressure considered is 143 psi. This 
pressure is the pressure required to cause a shell with a 
1.0-inch-long initial crack and a 0.040-inch wall thick- 
ness to fail due to an internal pressure load. The axial 
compression loads are increased from zero to the maxi- 
mum axial load that the shell can support. Typical re- 
sults are presented to illustrate the effects of crack length 
on the tearing response of a shell subjected to internal 
pressure loads, and on the prebuckling, buckling and 
postbuckling responses of a shell subjected to axial com- 
pression loads. Shells with both 0.020- and 0.040- inch 
wall thicknesses are considered for the analysis results, 
but only shells with 0.040-inch wall thickness are con- 
sidered for the experimental results. 

Internal Pressure Loads 

The effect of increasing the internal pressure in a 
shell on the total crack growth or crack extension is 
shown in figures 2a and 2b for shells with initial crack 
lengths of 1.0, 2.0, 3.0, and 4.0 inches, and with wall 
thicknesses of 0.040- and 0.020.-inches, respectively. 
These results indicate that the internal pressure in the 
shell can be increased, and the cracked shell will remain 
in equilibrium, up to a pressure at which yielding occurs 
at the crack tips and the opening angle at the crack tips 
reaches a critical value. At this pressure, the crack will 
start to grow. The initial growth of the crack is stable and 
the crack will not extend unless the pressure is increased. 
Eventually, unstable crack growth occurs. Unstable 
growth occurs when the slope of the curves in figures 2a 
and 2b becomes zero, which means that a small increase 
in pressure causes a very large increment in crack exten- 
sion. 

The effect of increasing the initial crack length on 
the internal pressure required to initiate stable tearing 
and on the pressure at which stable growth changes to 
unstable crack growth is summarized in figure 2c for the 
two shell wall thicknesses. The lower curve for each 
wall thickness represents the pressure required to initiate 
stable tearing as a function of normalized initial crack 
length, and the upper curve for each wall thickness rep- 
resents the pressure that causes the crack growth to be- 
come unstable. These results indicate that the internal 
pressure required to initiate stable tearing and to cause 
unstable crack growth decreases as the initial crack 
length increases. The results shown in figure 2 also indi- 
cate that the difference between the internal pressure re- 
quired to initiate stable tearing and the internal pressure 
that causes unstable tearing decreases as the initial crack 
length increases. For an initial crack length of 1.0 inch, 
stable tearing initiates in the 0.040-thick shell when the 
internal pressure is approximately 1 13 psi, and unstable 
crack growth occurs when the internal pressure is ap- 
proximately 143 psi. The difference between the internal 
pressure required to initiate stable tearing and the inter- 


nal pressure that causes unstable crack growth is approx- 
imately 30 psi for the 1.0-inch initial crack length in the 
0.040-inch-thick shell. The results for an initial crack 
length of 2.0 inches in the 0.040-inch-thick shell indicate 
that stable tearing initiates when the internal pressure is 
approximately 5 1 psi and unstable crack growth occurs 
when the internal pressure is approximately 79 psi. The 
difference between the internal pressure required to ini- 
tiate stable tearing and the internal pressure that causes 
unstable crack growth for this initial crack length is 28 
psi. The difference between the internal pressure re- 
quired to initiate stable tearing and the internal pressure 
that causes unstable crack growth for the 3.0- and 4.0- 
inch initial crack lengths in the 0.040-inch-thick shell is 
approximately 24 psi and 21 psi, respectively. The re- 
sults for the 0.020-inch-thick shell are similar, but with 
lower values of internal pressure. 

The results from the nonlinear analyses indicate 
that the internal pressure load induces large outward ra- 
dial displacements in the neighborhood of the crack. The 
response associated with these radial displacements is 
often referred to as “crack bulging” in the literature and 
is the cause of the larger crack opening displacements 
and crack-tip stress-intensity factors in a shell, compared 
to those for a flat sheet under otherwise identical condi- 
tions. 12 Furthermore, the extent of crack bulging is a 
function of the nondimensional crack curvature parame- 
ter, a/ Jrt. This observation is illustrated for the 0.040- 
inch-thick shell in figure 3, which shows the radial dis- 
placement response along the shell length at the circum- 
ferential location 0 = 0°, just before the initiation of 
stable crack growth, for each of the initial crack lengths. 
The radial displacement at the center of the crack nor- 
malized by the shell thickness, t , is w Q /t =1.00, 1.65, 
2.38 and 3.23 for the 1.0-, 2.0- 3.0- and 4.0-inch-long 
cracks, respectively. These displacements are greater 
than or equal to the shell wall thickness for all of the 
crack lengths considered, and represent large displace- 
ments in the context of nonlinear thin shell theory. The 
increase in crack bulging for the longer crack lengths is 
consistent with the observation that the pressure required 
to initiate stable tearing decreases with the increase in 
crack length. 

Experimental and analytical results for the 0.040- 
inch-thick shell are compared in figure 4 for shells with 
initial crack lengths of 2.0, 3.0, and 4.0 inches. As pre- 
viously mentioned, the experimental data obtained for 
the shell with a 4.0-inch-long initial crack were used to 
determine the critical CTOA angle to be used in predict- 
ing the response of the shells with 2.0- and 3.0-inch-long 
cracks. The analytical and experimental results correlate 
very well for the shells with the 2.0- and 3.0-inch-long 
initial crack lengths. 

Axial C ompression Loads 

The predicted load-shortening responses for a 
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0.040-inch-thick shell and a 0.020-inch-thick shell with 
initial crack lengths of 1.0, 2.0, 3.0 and 4.0 inches and 
subjected to axial compression are shown in figures 5a 
and 5b, respectively. An initial outward geometric im- 
perfection in the form of the lowest eigenmode was used 
in the nonlinear analyses to trigger local deformations in 
the vicinity of the crack. The applied compression and 
end-shortening values are normalized by the correspond- 
ing classical buckling values for a shell without a crack. 
For the 0.040-inch-thick shell and the 0.020-inch-thick 
shell with a 1.0-inch-long crack, the crack introduces an 
effective imperfection that causes general instability to 
occur at the loads indicated by the X in figures 5a and 5b. 
These shells cannot support additional compression load 
after buckling. For a shell with a longer crack, local 
buckling near the crack precedes shell collapse. The 
filled symbols in figures 5a and 5b identify the loads that 
correspond to initial local buckling near the crack for the 
0.040-inch-thick shell and the 0.020-inch-thick shell, re- 
spectively, with the 2.0-, 3.0-, and 4.0-inch-long initial 
cracks. Prior to buckling the radial displacement w Q at 
the center of the crack edges is nearly equal to zero. 
Once the critical load is reached, w 0 grows rapidly with 
increase in load. Initial local buckling is followed by a 
stable postbuckling response and the load can be further 
increased after local buckling has occurred near the crack 
edges. 

As the load is increased after initial local buckling 
has occurred, the 0.040-inch-thick shells with the 3.0- 
and 4.0-inch-long cracks, and the 0.020-inch-thick shells 
with the 2.0-, 3.0-, and 4.0-inch-long cracks experience 
a change in the local buckling mode. The initial post- 
buckling response of the shells after the mode change is 
unstable, and as a result, the axial load decreases after 
buckling occurs. The unstable transition region in the re- 
sponse predictions is indicated by the broken lines in fig- 
ures 5a and 5b. The unstable transition from the stable 
initial buckled configuration to the stable postbuckling 
configuration was determined by using the transient 
analysis capability in STAGS. The transient analysis 
was continued until the kinetic energy in the system was 
small. A time history of the kinetic energy during the 
transient analysis of the 0.040-inch-thick shell with a 
3.0-inch-long crack is shown in figure 6. The deformed 
shapes labeled A, B, C and D in figure 6 correspond to 
the points A, B, C and D on the load-end shortening and 
kinetic energy history curves, and show the development 
of the shell’s postbuckled response. Point A corresponds 
to the initial buckling deformation, points B, and C, cor- 
respond to solutions obtained during the transient analy- 
sis, and point D represents the stable postbuckled 
equilibrium state. Once a stable equilibrium state was 
determined from the transient analysis, the nonlinear 
static analysis was resumed to compute the stable post- 
buckling equilibrium response results shown in figure 
5a. The analysis was continued for the shells with the 


3.0- and 4.0-inch-long cracks for the 0.040-inch-thick 
shells, and for the shells with the 2.0-, 3.0-, and 4.0-inch- 
long cracks for the 0.020-inch-thick shells until compu- 
tational difficulties were encountered. The point of 
computational difficulty is assumed to correspond to 
shell collapse. Collapse of the 0.040-inch-thick and of 
the 0.020-inch-thick shells with 2.0-, 3.0- and 4.0-inch- 
long cracks is identified by the unfilled symbols in fig- 
ures 5a and 5b, respectively. 

The initial local buckling load predictions and qual- 
itative approximations for the shell collapse load for the 
0.020- and 0.040-inch-thick shells are summarized in 
figure 5c. The values of the normalized initial buckling 
loads for the 0.040-inch-thick shells are P/P cr = 0.88, 
0.54, 0.39, 0.32, and the qualitative approximate values 
of the normalized collapse loads are P/P cr = 0.88, 
0.59, 0.49, 0.47 for the 1.0-, 2.0-, 3.0- and 4.0-inch-long 
cracks, respectively. The values of the normalized initial 
buckling loads for the 0.020-inch-thick shells are 
P/P cr = 0.73, 0.41, 0.30, 0.24, and the qualitative ap- 
proximate values of the normalized collapse loads are 
P/P cr = 0.73, 0.47, 0.43, 0.41 for the 1.0-, 2.0-, 3.0- 
and 4.0-inch-long cracks, respectively. These results in- 
dicate that the magnitudes of the initial buckling loads 
and collapse loads for the shells decrease as the initial 
crack length increases. In addition, the difference be- 
tween the load at initial buckling and the load at collapse 
is smaller for the shorter crack lengths. The initial local 
buckling load results for the 0.020- and 0.040-inch-thick 
shells are represented very well by a characteristic curve 
that is based on the curvature parameter a/ Jrt. 

As the initial crack length increases, the radial dis- 
placements, w , near the crack increase. Radial displace- 
ment profiles along the shell length are compared in 
figure 7 for loads just before unstable buckling at the cir- 
cumferential location 8= 0°. The displacements are 
normalized by the shell wall thickness. The value of the 
radial displacement at the crack center, w 0 , normalized 
by the shell wall thickness, is w 0 /t = 0.2, 1.2, 2.1 and 
3.0 for the 1.0-, 2.0- 3.0- and 4.0-inch-long cracks, re- 
spectively. These prebuckling displacements are greater 
than the shell thickness for crack lengths greater than or 
equal to 2.0 inches, and represent large displacements in 
the context of nonlinear thin shell theory. 

Typical results of the elastic prebuckling, buckling, 
and postbuckling calculations are presented in figure 8 
for the 0.040-inch-thick shell with a 3.0-inch-long longi- 
tudinal crack. Resultant stress distribution contour plots 
obtained from the nonlinear analysis are shown on the 
corresponding deformed shapes in the figures. The dis- 
placement magnification factor for all images is five. In 
addition, on all contour plots, the darker regions of the 
figures represent higher values of the stress resultants, 
and the lighter regions represent lower values of the 
stress resultants. As described previously, the initial lo- 
cal buckling of the 0.040-inch- thick shells with the long- 
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er 3.0- and 4.0-inch~iong cracks is followed by a change 
in the local buckling mode and then a stable postbuckling 
response before the shell buckles into its general instabil- 
ity mode. A typical example of the initial local buckling 
displacement pattern for a 0.040-inch-thick shell with a 
3.0- inch-long crack is shown in figure 8a. The initial 
buckled deformation pattern, shown in figure 8a, is typi- 
cal of the shape of the initial buckled deformation pattern 
for the shells with the 2.0-, 3.0- and 4.0-inch-long cracks. 
The primary difference in the initial buckled deformation 
patterns for a shell with different crack lengths, with the 
exception of the 1.0-inch-long crack, is the magnitude of 
the deformations (see figure 7). The hoop and axial 
stress resultant distributions indicate that the crack has a 
significant and noticeable effect on the prebuckling 
stress distribution. The stable postbuckling deformation 
shape, after the local mode change, is shown in figure 8b. 
The deformation pattern in figure 8b indicates that the 
shell deforms into a local pattern with high circumferen- 
tial curvature that apparently stiffens the skin near the 
crack enough longitudinally to stabilize the shell and to 
increase the amount of axial compression load that can 
be supported by the shell after the local mode change oc- 
curs. 

The effects of material nonlinearities on the pre- 
dicted load-shortening response and deformation pat- 
terns or a 0.040-inch-thick shell with a 3.0-inch-long 
crack are shown in figure 9. Initial plastic yielding oc- 
curs for an applied load corresponding to P/P cr = 
0.43, which is approximately 90 percent of the buckling 
load at the local mode change. The load-shortening pre- 
dictions based upon an elastic analysis and an elastic- 
plastic analysis, and the deformation pattern prediction 
just prior to the local mode change, indicate that the shell 
response prior to the local mode change is adequately 
predicted by an elastic analysis. The initial postbuckled 
response predictions of the two analyses are, however, 
significantly different. The initial postbuckled deforma- 
tion pattern prediction from the elastic-plastic analyses is 
shown in figure 9c and the initial postbuckled deforma- 
tion pattern prediction from the elastic analysis is shown 
in figure 9d. Computational difficulties were encoun- 
tered in the elastic-plastic analysis at this point, prevent- 
ing prediction of the response further into the 
postbuckling load range. The postbuckled deformation 
pattern predicted by the elastic-plastic analysis is very 
similar to the deformation pattern labeled C in the elas- 
tic, transient analysis shown in figure 6. These results 
suggest that for the cylinder studied, yielding of the alu- 
minum will prevent the postbuckled deformation pattern 
predicted by the elastic analysis from developing. This 
suggestion is consistent with the experimentally ob- 
served behavior. The initial postbuckled deformation 
pattern for the 0.040-inch-thick cylinder with a 3.0-inch- 
long crack is very similar to the postbuckled deformation 
pattern predicted by the elastic-plastic analysis. 


The predicted prebuckling and initial buckling 
analysis results and the measured results for the 0.040- 
inch-thick shells with crack lengths of 2.0, 3.0, and 4.0 
inches are shown in figures 10a and 10b, respectively. 
Buckling of each of the shells is identified by the filled 
symbols in the figures, and the values of u/u cr and 
P/P cr at buckling are provided in parentheses. The 
analysis predictions were obtained using a material non- 
linear analysis. The measured loading-platen displace- 
ments were used to determine the amount of end 
rotations of the shell, and these data were used as input 
loading parameters for the nonlinear analysis. The re- 
sults shown in figures 10a and 10b indicate that the ana- 
lytical results slightly underestimate the experimental 
results for the initial buckling loads. The analytical un- 
derestimation of the initial buckling results could be a re- 
sult of uncertainties in the loading platen end rotations, 
and other differences in the analytical model and the as- 
tested specimens. An indication of the sensitivity of the 
predicted initial buckling load to variations in the input 
end rotations is provided in figure 10c. Figure 10c com- 
pares the measured load-shortening response for a shell 
with a 3.0-inch-long crack with analysis results that were 
obtained using two different values for the end rotation 
about one of the loading platen axes. These results indi- 
cate that the end rotations have some effect on the mag- 
nitude of the predicted buckling loads. 

Concluding Remarks 

The results of an analytical study of the effects of a 
longitudinal crack on the nonlinear response of thin, un- 
stiffened, aluminum cylindrical shells subjected to inter- 
nal pressure and axial compression loads are presented. 
The results indicate that the nonlinear interaction be- 
tween the in-plane stress resultants and the out-of-plane 
displacements near a crack in a thin shell can significant- 
ly affect the structural response of the shell. Large local 
stress and displacement gradients exist near a crack in a 
shell for all loading conditions considered in the study. 
The results indicate that the nonlinear response of a shell 
depends on the loading condition applied to the shell and 
the initial crack length. The magnitude of the internal 
pressure required to initiate stable tearing in a shell sub- 
jected to internal pressure decreases as the initial crack 
length increases. The magnitude of the internal pressure 
required to cause unstable crack growth in a shell also 
decreases as the initial crack length increases. The initial 
buckling load of a shell subjected to axial compression 
decreases as the initial crack length increases. Initial 
buckling causes general instability or collapse of a shell 
for shorter initial crack lengths. Initial buckling is a sta- 
ble local response mode for longer initial crack lengths. 
This stable local buckling response is followed by a sta- 
ble postbuckling response, which is followed by general 
or overall instability of the shell. 
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(b) Finite element model 
for pressure loads 



(c) Finite element model for axial 
, compression loads 


Figure 1. Shell geometry and finite element 
models. 
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(a) Crack extension for 0.040-inch-thick shells (b) Crack extension for 0.020-inch-thick shells 



a/ Jrt 

(c) Initial stable tearing and unstable crack growth 


Figure 2. Effect of increasing internal pressure on initial stable tearing and unstable crack 
growth for different initial crack lengths. 
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Figure 3. Radial displacement response along the shell length at the circumferential location 0 = 0°, just 
before the initiation of stable crack growth for 0.040-inch-thick shells. 



Total crack extension, Aa, in. 

Figure 4. Comparison of analytical and experimental total crack extension results for 0.040-inch-thick 
internally pressurized shells. 
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1.0 



(a) Load-shortening response for 
0.040-inch-thick shells 


l.Of 


X l-in.-long crack 

■ □ 2-in.-long crack 

• O 3-in.-long crack 

A A 4-in.-long crack 



(b) Load-shortening response for 
0.020-inch-thick shells 



(c) Initial buckling and collapse loads 

Figure 5. Effect of initial crack length on the response of cylindrical shells with longitudinal 
cracks and subjected to axial compression. 
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1.0 


Prebuckling, elastic 
and elastic-plastic analysis 



(a) Predicted load-shortening relation for linear-elastic 
response and material-nonlinear response 




(b) Predicted prebuckling deformation pattern for 
linear-elastic response and material-nonlinear 
response 



(c) Predicted initial postbuckling deformation < d ) Predicted initial postbuckling deformation 

pattern for material-nonlinear response pattern for linear-elastic response 


Figure 9. Comparison of linear-elastic and material-nonlinear responses of a 0.040-inch-thick shell with a 
3.0-inch-long crack and subjected to axial compression. 
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■ 2-in. -long crack • 3-in.-long crack A 4-in. -long crack 
Buckling values in parenthesis: (u/u cr , P/P cr ) 



(a) Predicted prebuckling load-shortening (b) Measured load-shortening relations 

relations 


■ Measured (0.46,0.45) 



(c) Effect of end rotations on the load-shortening results 
for a shell with a 3.0-inch-long crack 


Figure 10. Summary of predicted and experimental load-shortening relations for 0.040-inch thick shells 
subjected to axial compression. 
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